Excitonic energy transfer in light-harvesting complexes in purple bacteria 
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Two distinct approaches, the Frenkel-Dirac time-dependent variation and the Haken-Strobl model, 
are adopted to study energy transfer dynamics in single-ring and double-ring light-harvesting sys- 
tems in purple bacteria. It is found that inclusion of long-range dipolar interactions in the two 
methods results in significant increases in intra- or inter-ring exciton transfer efficiency. The depen- 
dence of exciton transfer efficiency on trapping positions on single rings of LH2 (B850) and LH1 is 
similar to that in toy models with nearest-neighbor coupling only. However, owing to the symmetry 
breaking caused by the dimerization of BChls and dipolar couplings, such dependence has been 
largely suppressed. In the studies of coupled-ring systems, both methods reveal interesting role 
of dipolar interaction in increasing energy transfer efficiency by introducing multiple intra/inter- 
ring transfer paths. Importantly, the time scale ( 4ps) of inter-ring exciton transfer obtained from 
polaron dynamics is in good agreement with previous studies. In a double-ring LH2 system, dipole- 
induced symmetry breaking leads to global minima and local minima of the average trapping time 
when there is a finite value of non-zero dephasing rate, suggesting that environment plays a role 
in preserving quantum coherent energy transfer. In contrast, dephasing comes into play only when 
the perfect cylindrical symmetry in the hypothetic system is broken. This study has revealed that 
dipolar interaction between chromophores may play an important part in the high energy transfer 
efficiency in the LH2 system and many other natural photosynthetic systems. 

PACS numbers: 



I. INTRODUCTION 

An essential photochemical process in plants and pho- 
tosynthetic bacteria is light-harvesting, in which solar 
photons are first absorbed by antenna complexes (also 
known as the antenna protein), and the resulting exci- 
tonic energy is transferred with a high efficiency to spe- 
cialized pigment-protein complexes (PPC) or the reac- 
tion centers (RC). In a photosynthetic system of bacteri- 
ochlorophylls aggregates embedded in protein scaffolds, it 
has been suggested that the protein environment plays a 
critical role in protecting excitonic coherence, and conse- 
quently, facilitating coherent energy transfer. The entire 
process of coherent energy transfer takes place in a very 
short period of time, from a few hundred femtoseconds to 
a few picoseconds. Recent studies have shown that the 
efficiency of energy transfer from antenna pigments to 
the reaction centers is remarkably high (e.g., >95%)i~— . 
Such a high efficiency of energy transfer has inspired the 
designs of efficient artificial systems to convert solar en- 
ergy into other forms. However, an acute lack of reli- 
able simulation techniques and theoretical tools for un- 
derstanding the exact mechanisms of long time quantum 
coherence in such systems hinders the realistic applica- 
tion of the highly efficient energy transfer in designing 
artificial systems. Currently, the Forster resonance en- 
ergy transfer theory (FRET) based on weak electronic 
coupling approximation and the Redfield equation based 



on Markovian approximation are the two commonly ap- 
plied methods to deal with the energy transfer processes 
in these complexes. The classical Forster theory treats 
electronic coupling perturbatively and also assumes inco- 
herent hopping of exciton between single donor and ac- 
ceptor induced by point dipole-dipole interaction of tran- 
sition dipoles 5 of the chromophores. Assumption of point 
dipole is inadequate in describing pigments' aggregation 
which results in the breakdown of the approximation 6 . 
To overcome the drawback of the point-dipole approx- 
imation, Jang et all ~— developed a muti-chromophorc 
version of FRET theory which can be applied to cope 
with the transfer pathway interference that is not ac- 
counted for in the previous FRET treatments. However, 
when the electronic coupling is strong but the system- 
bath coupling is weak, it is necessary to consider relax- 
ation between delocalized exciton states. In this limit, 
the dynamics of exciton energy transfer (EET) can be 
described by the coupled Redfield equations in the exci- 
ton basi a 12 ' 13 ' 15 ' 16 . Redfield equation, on the other hand, 
is applicable in the weak exciton-phonon coupling limit 
where the problem can be treated perturbatively. More 
importantly, in the Redfield equation-^, the bath degree 
of freedom has been projected out to obtain a reduced 
master equation for the density matrix evolution under 
the Markovian approximation. However, for a function- 
ing photosynthetic system, the exciton-phonon coupling 
is not weak and also Markovian approximation is a poor 
one in describing the realistic bath behavior since the 
typical phonon relaxation time scale is quite long. 
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In purple bacteria, the antenna complexes have 
highly symmetric multi-ring structure s 20 ' 21 , which are 
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the subject of numerous studies regarding the mech- 
anism underlying the highly efficient energy transfer 
processe d 3 ' 4 ' 23 " — . Photons are absorbed mainly by 
the peripheral antennae complex LH2 of purple bac- 
teria, where the photo-generated excitons are sub- 
sequently transferred to the reaction center within 
the harvesting complexe LH1 in a few picoseconds—. 
The LH2 complex is an eight-unit circular aggregate 
built from a/3-heterodimers forming C8 symmetr y 20 ' 21 . 
Each unit contains a pair of a and ft apoproteins, 
three bacteriochlorophylls-a(BChls-a) molecules and a 
carotenoid, and the BChl-a molecules form two rings 
named according to their corresponding absorption max- 
ima at 800nm and 850nm as the B800 and B850 rings, 
respectively. The B850 ring consists of 16 tightly po- 
sitioned BChls-a, with the Mg-Mg distance about 9.36 
A for the la — 1(3 dimer, and about 8.78 A for the 
la — 1/3 dimer~. There are many studies on the 
mechanism behind these highly efficient photosynthetic 
processe s 3 ' 4 ' 23- — . In Refi^S, optimal energy transfer con- 
ditions in linear and nonlinear systems were analyzed, 
and a similar approach was recently applied to self- 
assembling light-harvesting arrays^ to study optimal ef- 
ficiency. For the antenna complexes of purple bacteria, 
energy transfer dynamics in LH1 and LH2 rings is inad- 
equately understood, and questions remain on regarding 
how dipolar interaction, symmetry of the ring and inter- 
ring coupling etc., affect the energy transfer efficiency. 

Inspired by previous studies, in this work we study in 
detail the effect of intra-ring dipolar interactions, ring 
geometry and symmetry, and more importantly, inter- 
ring distances in determining EET. As a starting point 
of understanding EET events, two approaches have been 
adopted to investigate various light-harvesting complexes 
in purple bacteria. For the low-temperature exciton dy- 
namics, we use the Frenkel-Dirac time-dependent vari- 
ational method to simulate the quantum dynamics of 
the Holstein polaron via the Davydov Dl Ansatz, while 
for the high-temperature dynamics and transfer efficiency 
calculations, the Haken-Strobl model 1 ^, equivalent to the 
Redfield equation at high temperatures, has been em- 
ployed. 

This paper is organized as follows. In section II, the 
methodology of aiming for establishing possible optimal 
conditions and understanding of dynamics for energy 
transfer in the ring systems has been proposed. The ef- 
fect of various parameters including the effect of dynam- 
ics disorder due to coupling to phonons and dephasing 
due to inclusion of environment in determining optimal 
energy transfer efficiency has been discussed in details in 
section III. Finally the conclusions are drawn in section 
IV. 



II. METHODOLOGY 

Two approaches are adopted to study the energy trans- 
fer process in the light-harvesting apparatus of the pur- 



ple bacteria, the Frenkel-Dirac time-dependent varia- 
tional method making use of the Holstein model which 
can be generally regarded as low-temperature cases, and 
the Haken-Strobl method applicable to high-temperature 
scenarios. 

A. Polaron Dynamics in LH2 Systems 

The Holstein Hamiltonian for the exciton-phonon sys- 
tem reads a o 30 ' 31 

H = H cx + H p h + i/cx-ph (1) 

with 

H ex = - J ^ QnQWl + ttn-l); ( 2 ) 



H ph = ^2^qblbq, (3) 
9 

H cx - ph = -~7=J2J2 9qu q ala n (b q e iqn + H.c.) , (4) 

* n q 

here a) n (a n ) is the creation (annihilation) operator for 
an exciton at the n th site, and J is the exciton transfer 
integral. b q (b q ) is the creation (annihilation) operator of 
a phonon with momentum q and frequency uj q , and the 

Planck's constant is set as H = 1. i? cx _ p h is the linear, 
diagonal exciton-phonon coupling Hamiltonian with g q 
as the coupling constant, with J^qdfaq = ScJo- Here S 
is the well-known Huang-Rhys factor, and u>o is the char- 
acteristic phonon frequency, which is set to unity in this 
paper for simplicity. It is noted that in this method, 
phonon dispersions in the Hamiltonian can have vari- 
ous forms to describe many phonon branches of different 
origins. For example, a linear phonon dispersion with 
u q = lj [1 + W(2\q\/ir — 1)] may be adopted, where W is 
a constant between and 1, and the band width of the 
phonon frequency is 2Wujo. Furthermore, other forms 
of exciton-phonon coupling can also be included, which 
implies the capability of current method in imitating the 
spectral density of phonon obtained with other methods 
such as molecular dynamics 32 . 

For multiple-ring system, by assuming only exciton 
transfer integral and dipolar interaction between rings, 
the original Holstein Hamiltonian can be modified as: 

^cx = — y ] y ] J nm&n^ Q'rn-i (5) 
rir2 nm 

where n = l...N x i ag and = l..-/V r j n g indicates any one 
pair of rings (with total number of N r i ng rings) within the 
2-dimensional lattice if ri ^ r2, where the summation is 
over all rings. However if n = r-z, Eq. ([5]) becomes: 
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where r runs from 1 to N r i ng rings. For each ring, the 
exciton Hamiltonian takes the form of Frenkel exciton 
Hamiltonian with J nm given by 36 : 
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where E\ and £2 are the on-site excitation energies of an 
individual BChl-a, J\ and J 2 are the transfer integral 
between nearest neighbors, N equals 8 as the system is 
of C8 symmetry, and matrix elements Wij are the dipo- 
lar coupling for the non-nearest neighbors. The dipolar 
coupling between site i and site j takes the form 



W itj = C 



<1- <I , 3(di • Tij)(dj ■ Tij) 



(8) 



where C is the proportionality constant, is the vec- 
tor connecting the ith and jth monomers, and are the 
unit vectors of the ith BChl-a. In our calculations in 
this work, the following parameters^ 3 - are adopted: J\ = 
594cm- 1 , J 2 = 491cm- 1 and C = 640725A 3 c m - 1 . When 
the J n m matrix is applied to the polaron dynamics cal- 
culations, it must be scaled by the characteristic phonon 
frequency ujo- The Mg-Mg distance between neighboring 
B850 BChls is 9.2 A within the a/3-heterodimer and 9.25 
A between neighboring heterodimers. For simplicity, the 
a and /3 proteins are assumed to be evenly distributed 
with a distance of 9.25 A. 

For multiple-ring systems, the phonon and exciton- 
phonon interaction Hamiltonian can be described as: 



EE< 



(9) 



H, 



cx— ph 



E E E 9 r q ^a r n (e lqn b r q + H.c 



( 10 ) 

where r is summed over 1 to -/V r ; ng . b q * {b q ) is the cre- 
ation (annihilation) operator of a phonon with momen- 
tum q and frequency uf q in the r ring. Similar to single- 
ring case, we can obtain the Huang-Rhys factor 5* from 
^2q(9q) 2uj q = Sujq. It is clearly shown in the formalism 
that the phonons in different LH2 rings are completely 
independent of each other in the current treatment. 

The family of the time-dependent Davydov Ansatze, 
i.e., the Di, D and D 2 trial states, can be written in a 
general form as 



|*n(t)> = X; a »(*)^(*)|0>-x|0) P h, 



(11) 



where a n (t) are the variational parameters representing 
exciton amplitudes, and Un(t) is the Glauber coherent 



operator 



^)=exp{£[A„,<r(t)&t-H.c.]}. 



(12) 



For the Di Ansatz, A„ i9 (i) are NxN variational param- 
eters representing phonon displacements. The D 2 and D 
Ansatze are two simplified cases of the Di Ansatz. In the 
D 2 Ansatz, X n ,q(t) are replaced by N independent vari- 
ational parameters f3 q (t), i.e., X n ,q(t) — P q (t). While in 
the D Ansatz, X n , q {t) are replaced by 2N — 1 variational 
parameters X (t), /3 q ^o(t) and 7g/o(*)> where X niq=0 (t) = 
A (t), and X n , q (t) = (3 q {t) + e-^ n lq (t), q^O. 

The time evolution of the photo-excited state in a 
one-dimensional molecular aggregate follows the time- 
dependent Schrodinger equation. There are several ap- 
proaches to solve the time-dependent Schrodinger equa- 
tion. In the Hilbert space, for example, the time- 
dependent wave function |$(f)) for the Hamiltonian H is 
parameterized by a set of time-dependent variables a m (t) 
(to = 1, ...,M): 



|S(f)> = \{a m (t)}). 



(13) 



Assuming that \&{t)) satisfies the time-dependent 
Schrodinger equation, one has 



(14) 



Explicitly putting in the Hamiltonian H of Eq. (fTJ and 
writing d\$(t))/dt in terms of a m (t) and their time- 
derivatives a m (t) (m = 1, M), one obtains 

- H)Mt)) = \{a m (t)}, \{a m (t)}) = 0- (15) 

Projecting Eq. (TIB"]) onto M different states l^m) (to — 
1, M), one obtains M equations of motion for the pa- 
rameters a m (t): 



(* m \{a m (t)},{a m (t)}) = 0. 



(16) 



The approach we adopt in this work is the Lagrangian 
formalism of the Dirac-Frenkel time-dependent varia- 
tional method^, a powerful technique to obtain approxi- 
mate dynamics of many-body quantum systems for which 
exact solutions often elude researchers. We formulate the 
Lagrangian L as follows 



L = <*(*)! 



ih d 
~2~dt 



Am))- 



(17) 



From this Lagrangian, equations of motion for the M 
functions of time, parameters a m {t), and their time- 
derivatives a m (t) (to = 1, M), can be obtained by 



d dL 
dt da* 



dL 

do.* 



0. 



(18) 
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Furthermore, to better elucidate the correlation between 
the exciton and the phonons, the phonon displacement 
£ n (t) is introduced and defined as 



£ n (i) = (tt D (i) 



Si 



*d(*) 



(19) 



where 61 



iqnij 
u q . 



The reader is referred to 



Appendix A for discussions on the precision of the Davy- 
dov trial states (defined as A(t)) and Appendix B for 
detailed derivations of equations of motion for polaron 
dynamics in a multiple-ring system. 



B. The Haken-Strobl model 

The second approach applied to the LH1/LH2 com- 
plexes in this work is the Haken-Strobl model, which has 
been previously used to describe the exciton dynamics 
in We start with the Liouville equation for the 

exciton density matrix p 

p(t) = -C P (t) 

[^sys ~t~ ^dissip 4~ ^dccay ^trapjp(^)- (^0) 

Here £ sys is the free evolution operation of the pure ex- 
citon system, and can be written as C sys = i[H,p]/h, 
where [H] nm = (1 - 5 nm )J nm + S nm e n with J nm the ex- 
citonic coupling between site n and site to, and e the 
site energy. £dissip represents the dephasing and popu- 
lation effects within the exciton manifold. For the con- 
venience of analysis, the coupling to the environment is 
given by the standard Bloch-Redfield equation, which can 
be written as [£<iissip]nm = (1 - S nm T nm ), where T is 
the pure dephasing rate. £dccay represents the decay of 
the exciton to the ground state and can be expressed as 
[£dccay]nm = (kd,n + k^ m )/2 with k^ n the decay rate on 
site n. £trap represents the trapping rate, which is criti- 
cal to the energy absorbing process, and can be described 
as [£trap]nm = (fa,n + fet.m)/ 2 with fa,n tne trapping rate 
on site n. 

In fact, among £dccay and Arap, which represent two 
possible channels for irreversible exciton energy loss, the 
former is ineffective in the energy transfer. Efficiency of 
energy transfer can be gauged by the quantum yield q, 
namely the trapping probability^!, 



X/n kt,n T n + X/n kd,n T n 



(21) 



where the mean residence time r„ can be expressed as 
r n = Jq 00 p n (t)dt, and the population p n — p nn is the di- 
agonal element of the density matrix. In photosynthetic 
systems, k^ is on the order of ps and fc^ 1 is on the order 
of ns, so the trapping rate is much larger than decay rate, 
and the quantum yield is close to unity. The quantum 
yield can then be approximated as q « (1 + fcd(i)) -1 , 
where (t) = r n is the mean first passage time to the 
trap state without the presence of the constant decay, 



i.e., the average trapping time. Quantum yield and trap- 
ping time have been extensively studied in the context of 
molecule photon statistics^ and experimental investiga- 
tion was also carried out in a photosynthetic system 3 ^. A 
shorter trapping time corresponds to a higher efficiency 
of exciton transfer. 



III. RESULTS AND DISCUSSIONS 
A. Polaron Dynamics in LH2 Systems 

With the aid of Frenkel-Dirac time-dependent varia- 
tional method, a hierarchy of variational wave functions 
has been previously formulated to faithfully describe the 
dynamics of the Holstein polaro n 37 i 38 . The accuracy of 
the Davydov Ansatze has been discussed in length in 
Ref i 37 ' 38 , and it has been shown that for diagonal exciton- 
phonon coupling, the Dl Ansatz has the highest precision 
in the hierarchy of the Davydov trial states when applied 
to the solution of Holstein Hamiltonian. The increase in 
the precision of the Dl Ansatz is achieved mainly by im- 
proved sophistication of the phonon wave functio n 37 ' 38 . 
Such increased sophistication plays a dramatic role espe- 
cially in the weak and moderate coupling regime, where 
the plane- wave-like phonon components render impracti- 
cal to use a simple linear superposition of phonon coher- 
ent states. The vastly increased number of variational 
parameters for the phonon wave function ultimately re- 
sults in higher precision of the Dl Ansatz. 

The approach is first to be applied to the exciton po- 
laron evolving in an isolated LH2 ring. To assure the 
accuracy of our approach, the Dl Ansatz has been cho- 
sen to deal with the dynamics in this paper. For a case 
of wide phonon bandwidth (W=0.8) and a moderate 
exciton-phonon coupling strength (S=0.5), time evolu- 
tion of exciton probability |a„(i)| 2 and phonon displace- 
ment (t) for the case with only nearest neighbor trans- 
fer integral have been displayed in Figs. [IJa) and (b). 
The polaron state is initialized as follows: a n (0) = 5 n ^ 
and A nj g(0) = 0. The corresponding dynamics with real- 
istic LH2 transfer integral matrix (c/., Eq. ((7|)) has been 
given in Figs. [2(c) and (d). Time-dependent exciton 
and phonon energies as well as the magnitude of devi- 
ation vector A(t) for both cases are plotted in Fig. [lje). 
For all calculations, a characteristic phonon frequency of 
luq = 1670cm -1 and a Huang-Rhys factor of S=0.5 are 
chosen^. The phonon frequency obtained from molecu- 
lar dynamics simulation is related to the stretching mode 
in either a C=0 bond or a methine bridge. 

It is important to compare calculated absorption spec- 
tra with measured ones in order to further validate our 
approach and lend support to our simulation of excita- 
tion transfer dynamics in a photosynthetic complex. Lin- 
ear absorption spectra are calculated, and compared to 
measured ones in Fig. [5J Detailed derivations for ab- 
sorption spectral calculations can be found in Refi 37 ' 38 
Measured data for B850 complexes are obtained from 
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FIG. 1: Real space dynamics of (a) exciton probability 
|an(t)| 2 and (b) phonon displacement £n(i) for a single-ring 
with only nearest neighbor transfer integral J = 0.3557, a 
moderate exciton-phonon coupling strength S = 0.5, and a 
wide bandwidth W = 0.8 of phonon dispersion. Correspond- 
ing real space dynamics are displayed for (c) exciton proba- 
bility |a n (i)| 2 and (d) phonon displacement for a single 
realistic LH2 ring transfer integral matrix. The Dl Ansatz is 
used for the simulation. The time evolution of energy com- 
ponents for exciton, phonon, exciton-phonon interaction and 
deviation magnitude A(f) for both cases has also been given 
in (e). Note that the value of the nearest neighbor transfer 
integral is taken from the largest element (scaled by the char- 
acteristic phonon frequency u>o) of the realistic LH2 transfer 
integral matrix. 



Ref4S Good agreement is found between theory and ex- 
periment, pointing to the applicability of our model to 
the B850 ring for the parameter set chosen. It is espe- 
cially worth noting that the long tail shown in the mea- 
sured spectrum is reproduced by our approach. 

Owing to the aforementioned higher accuracy of the Dl 
Ansatz as recorded in Refi 37 ' 38 in comparison with other 



FIG. 2: Comparison of calculated linear absorption spectra 
and measured ones. A single-ring with realistic LH2 ring 
(B850) transfer integral matrix is adopted. Parameters used 
are S = 0.5, u)q — 1670cm -1 , W = 0.8 and a damping factor 
of 0.12. 



trial states, as shown in Fig. fflje), it is therefore an ideal 
candidate for the study of excitation dynamics of the LH2 
system. The real LH2 system contains pigments that 
interact via long-range dipolar couplings as given by the 
pure exciton Hamiltonian with transfer integral matrix 
given in Eq. (|7j) . By introducing dipole-dipole interaction 
into the Hamiltonian, the accuracy for the Dl Ansatz 
has to be rechecked and compared to that for a model 
system with nearest neighbor coupling only. A first step 
toward the accuracy checkup is the comparison of polaron 
dynamics in Figs.QJa), (b), (c) and (d), where dramatic 
changes can be found as a result of dipolar interactions. 

In addition to the case with only nearest neighbor J 
shown in Figs. QJa) and (b), many new features have 
emerged as shown in Figs. HJc) and (d). One major 
finding is that the inclusion of realistic dipolar interac- 
tions in the LH2 system breaks the symmetry of the ring 
structure, resulting in asymmetric and faster exciton dy- 
namics compared to the model system with only nearest 
neighbor J. Such an effect can be related to an emer- 
gence of additional energy transfer paths when dipolar 
interactions are included. Although the dynamics here 
does not contain the coupling to the external heat bath, 
the essential physics is captured to a great extent. It 
is worth noting that in Figs. Hfb) andQJd), time evolu- 
tion of phonon displacement components, i.e., £, n {t) with 
"V" -shaped feature survived as the exciton Hamiltonian 
has been completely changed from nearest neighbor J to 
a realistic J nm - Which can be explained as follows: the 
phonon dynamics of the system is not only determined 
by the coupling between exciton and phonon, but also 
has its own propagation pattern which is determined by 
the linear dispersion relationship given earlier, thus once 
the excitation on site 8 starts to propagate as shown in 
the exciton dynamics, the strong lattice distortion in the 
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vicinity of the excitation will also propagate. Moreover, 
the effect of the exciton on phonon dynamics is apparent 
by comparing Figs, [lja) and (b), where an exact cor- 
respondence between the exciton dynamics and phonon 
can be found. However, the value of exciton-phonon cou- 
pling strength S tells us the coupling between exciton and 
phonon is moderate, thus the exciton generated on site 8 
can easily propagate through larger distance similar to a 
bare exciton. 

More interestingly, the revival and localization phe- 
nomenon shown in Fig. fflja) is nearly absent in Fig. [2c ) 
as a result of broken symmetry as well as multiple trans- 
fer channels due to the presence of dipolar interaction 
within the ring, although there are certain possible po- 
sitions such as on site 8 at around 6 phonon periods ap- 
pears to have possible revival of exciton amplitude sim- 
ilar to those shown in Fig. [Ua) on site 8 at around 10 
phonon periods, lending further support to the speed up 
of exciton dynamics (in terms of derealization) due to 
the inclusion of dipolar interaction. From above discus- 
sions, by including long-range excitonic coupling other 
than nearest neighbor (which is true for many realistic 
cases such as the LH2 system investigated in this paper), 
the exciton can more efficiently and evenly spread over 
the entire ring in the presence of exciton-phonon cou- 
pling. In the other words, nature has designed a unique 
way for the light harvesting complexes such as LH2 in 
maximizing their excitonic energy transfer efficiency in 
spite of exciton-phonon coupling responsible for localiza- 
tion and dephasing of the excitons. It appears the sym- 
metry breaking of the energy transfer pathway has more 
positive role in LH2 complexes. 

Furthermore, It is also important to note in Fig. [TJe) 
that for the chosen parameter sets, the Dl Ansatz is ca- 
pable to describe the realistic LH2 system with great 
precision. Comparing to the maximum absolute value 
of phonon or exciton-phonon interaction energy shown 
in the figure, the magnitude of deviation vector A(t) is 
less then 30 percent of the average system energy (here 
we use phonon energy as the reference). Inclusion of 
complex dipolar interaction does not further increase the 
error of the Dl Ansatz in describing the system dynam- 
ics. Although the deviation magnitude does not change 
much, the dynamics of exciton and phonon deviate dra- 
matically if one compares Figs. [1] (a), (b) andQ](c), (d). 

Using the Dl Ansatz, we also examined exciton- 
phonon dynamics of a realistic two-ring system with long- 
range dipolar interactions. Results from the dynamics 
calculations are shown in Fig. [3] and Fig. 0] for two inter- 
ring distances (15 and 20 A, respectively). The choice 
of the distances for the two-ring system is in accordance 
with the one shown in Refill, which suggests the Mg-Mg 
distance of around 24 A for the closest chromophores be- 
tween two adjacent LH2 rings. Considering the size of 
the chromophores, it is beneficial to slightly reduce the 
distance to consider the exciton transfer from the tip of 
one molecule to the other, thus two different distances, 
i.e., 15 and 20 A are used in this study. To have bet- 
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FIG. 3: Real space dynamics of (a) exciton probability 
|Qn(t)| 2 and (b) phonon displacement £ n (i) in ring-1 for a 
two-LH2-ring system with realistic LH2 ring transfer integral 
matrix, a moderate exciton-phonon coupling strength S = 0.5 
and a wide bandwidth W = 0.8 of phonon dispersion. The Dl 
Ansatz is used for the simulation. Corresponding real space 
dynamics in ring-2 are displayed for (c) exciton probability 
|o!n(t)| 2 and (d) phonon displacement £n(i). The time evo- 
lution of energy components from exciton, phonon, exciton- 
phonon interaction, deviation magnitude A(f) has also been 
given in (e). The distance between two rings is set to 15 A, 
with interaction between two rings given by Eq. (|5J. (f) The 
schematics of two connected LH2 complexes, r is the distance 
between the tips of the rings. 



7 





[i 0.0ns 0.01 oojs 0.03 0.025 n.oj -0.0! -0.005 0,005 o.ni 0.015 0.02 

fc 1 1 1 1 1 1 1 




5 5 10 15 20 25 



FIG. 4: Real space dynamics of (a) exciton probability 
|a n (i)| 2 and (b) phonon displacement £ n (t) in ring-1 for a 
two-LH2-ring system with realistic LH2 ring transfer integral 
matrix, a moderate exciton-phonon coupling strength S = 0.5 
and a wide bandwidth W = 0.8 of phonon dispersion. The Dl 
Ansatz is used for the simulation. Corresponding real space 
dynamics in ring-2 are displayed for (c) exciton probability 
|a„(t)| 2 and (d) phonon displacement £, n (t). The time evo- 
lution of energy components from exciton, phonon, exciton- 
phonon interaction, deviation magnitude A(t) has also been 
given in (e). The distance between two rings is set to 20 A, 
with interaction between two rings given by Eq. ©. 



ter understanding of the roles of intra- and inter-ring 
excitonic transfer, the two-ring system is initialized as 
follows: o£(0) = 5r,i^ri,8 and A£ g (0) = 0, i.e., a sin- 
gle excitation at site 8 of ring-1, after a subsequential 
spreading within ring-1, the exciton population is trans- 
ferred to site 1, which is the closet point to site 10 (and 
nearby sites) of ring-2, at about 3 phonon periods from 
t=0. This exciton transfer process is clearly revealed in 
Figs. Hie) and He). 

It is always helpful to compare the value of A(t) rela- 
tive with system energies as the first check. As shown in 
Figs.|3fe) andHJe), the accuracy of the Dl Ansatz applied 



to the two-ring system with different inter-ring distances 
does not alter much as compared to the one-ring sys- 
tem since the magnitude A(t) remains close to that for 
the one-ring case. The value of A(t) also becomes fairly 
stable after initial 5 phonon periods, therefore confirms 
that the method applied to study the multiple-ring ex- 
citon dynamics using time-dependent variational Ansatz 
is valid. 

At a first glimpse, one may spot few differences be- 
tween Fig. [3] and Fig. U during the first 10 phonon peri- 
ods. One can also find resemblance of the exciton prob- 
abilities and phonon displacement between the two-ring 
system and the one-ring system in Fig. [T] if comparison 
is made during the same period of time. This is due to 
the fact that a very small amount of exciton population is 
transferred to the second ring of the two-ring system dur- 
ing this period as clearly revealed in Figs. |3Jc) andSJc), 
thus the dynamics of both the exciton and phonons are 
essentially those of the one-ring case. However, after 10 
phonon periods, as the transfer of exciton population in- 
creases, |a„(i)| 2 and £ n (i) become quite different for the 
two cases of different inter-ring distances. 

Interestingly, as one may observe from Figs. 0D[c) 
and[4jc), many bright spots appear in the exciton prob- 
ability contour of ring-2, indicating possible pathway in- 
terference within the ring. The origins of the bright spots 
are explained as follows. Among the nearest neighbors of 
ring-1 [cf. Fig.[3](f)], sites 8 to 10 of the second ring show 
more pronounced changes in exciton probability, if one 
compares Figs. (He) andHJc). It is also true that inter- 
actions between sites 1, 2 and 16 of ring-1 and sites 8, 9 
and 10 of ring-2 dominate inter-ring transfer according to 
Eq. ([7]), resulting in multiple inter- ring exciton transfer 
pathways. 

Due to strong interactions between site 1 of ring-1 and 
site 9 of ring-2 owing to their proximity, it is somewhat 
counterintuitive to observe the relatively smaller exciton 
probability on site 9 of ring-2. The explanation lies in the 
dimerization of LH2 rings which results in opposite signs 
of dipole moments of site 1 of ring-1 and site 9 of ring- 
2, as shown in Fig. I^f). Besides dimerization, it is also 
helpful to consider the intra-ring interference in ring-2 be- 
tween exciton wave packets that travel through multiple 
paths from ring-1 to ring-2. Figs. (He) andSJc) suggest 
that sites 8 and 10 of ring-2 (especially, site 10) can be 
treated as input sites strongly connected to ring-1. Ap- 
pearance of bright spots on sites 8 and 10 can be simply 
understood as a combined effect of exciton transfer from 
ring-1 to ring-2 and constructive interference of exciton 
wave packets in ring-2. However, on locations far from 
site 10, such as sites 2, 5, 13 and 15 of the second ring, 
the effect of constructive interference is more prominent. 
As sites 2 and 16 of ring-2 are opposite to the aforemen- 
tioned input sites, it is not surprising to observe bright 
spots on site 2 and 16 as a result of constructive interfer- 
ence, as shown in Figs.[3][c) and0|c). Effects of multiple 
pathways for inter-ring exciton transfer may play a role 
in the appearance of additional bright spots on ring-2. 
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FIG. 5: Exciton population dynamics (Y] n |a^(t)| 2 ) for the 
two different separations shown in Figs. [2] and O where a 
long time dynamics of the exciton population has also been 
given in the inset for the ring system with a distance of 15 A. 
Exciton population dynamics for the inter-ring distance of 15 
A with on-site energy static disorder has also been given in 
the figure. 



Even more interestingly, one can also calculate the ex- 
citon population of a given ring, i.e., J2 n l a nW| 2 m rm g- 
r (with r the ring index) as a function of time to shed 
light on the overall exciton dynamics of the two-ring sys- 
tem with inter-ring distance of 15 and 20 A. Calculation 
results are displayed in Fig. [5l where exciton populations 
corresponding to Figs. [3] and 2] are plotted as a func- 
tion of time. A supplementary calculation for a much 
longer time is also given in the inset of Fig. [5] for inter- 
ring distance of 15 A. Short-time results suggest a nearly 
stepwise population transfer from ring-1 to ring-2. The 
stepwise features of population profiles in Fig. [5] coincide 
with the increases of exciton probabilities on sites 8 to 10 
in the second ring as shown in Fig.[3fc) and@|c), lending 
further support to our previous discussion on inter-ring 
exciton transfer. Long-time calculation for the case with 
an inter-ring distance of 15 A reveals that it takes about 
200 phonon periods to have complete transfer of the ex- 
citon population, as shown in the inset of Fig. [5] As one 
phonon period is equal to about 20/s, the estimated to- 
tal transfer time of 4ps is in good agreement with the 
findings in Refj2£. 

Effects of static disorder on polaron dynamics have also 
been examined in order to test the applicability of the 
current approach to photosynthetic systems. Simulation 
results with Gaussian-distributed static disorder are in- 
cluded in Fig. [5j The on-site energy average is set to 
zero, and a standard deviation of a = O.lwo is adopted 
in accordance with Ref.— . For simplicity, only on-site en- 
ergy disorder is included although it will not be difficult 
to examine off-diagonal disorder. Over 100 realizations 
have been simulated in the calculation. Fig. [S] provides a 



comparison on the exciton population transfer with and 
without the energetic disorder. While some population 
fluctuations have been smoothened out owing to disorder, 
the overall energy transfer efficiency is not significantly 
affected by the addition of static disorder. Furthermore, 
in the absence of any exact or more reliable dynamics 
studies, our approach here provide useful guide on the 
dynamics of excitonic energy transfer in LHI-LHII sys- 
tems. 



B. Exciton Population Dynamics of N-site 
Symmetric Simple Ring Systems and LH2 
Complexes 

Dynamics study in the previous section reveals the ef- 
fect of dipolar interactions on exciton transfer efficiency 
at low temperatures, but it is also interesting to look at 
the opposite limit of high temperatures. As a starting 
point, it is preferable to use simple ring systems to ana- 
lyze what may affect the energy transfer efficiency at high 
temperatures. Here we first study the energy transfer of 
symmetric ring systems of three, four, six, seven, eight 
and sixteen pigments. Efficiency of the energy transfer 
with respect to various trapping positions will be exam- 
ined, and for an iV-site ring system, analytical solutions 
will be obtained using the stationary condition on the 
off-diagonal density matrix element p nm . 

For a three-site ring system, we assume that all the site 
energies are the same, thus the energy difference is zero. 
The coupling between neighboring sites is assumed to be 
constant such that Jn = J 2 i = J13 = J31 = J23 = J32 = 
J. Furthermore, the coupling with the environment (de- 
phasing) is Ti2 = r 2 i = Ti 3 = r 3 i = r 23 = r 32 = r. 

If the system is excited at site 1 and the excitation is 
trapped at the same location, the average trapping time 
is given by (t) = 3fc t _1 , independent of V and J. If a pho- 
ton is absorbed at site 1, and the excitation is trapped 
at site 2, the average trapping time takes the form 

3 1 1 2T + k t 

{t) -k t + f + 2TTk t + ^ (22) 

which depends on T and J. It is greater than that of the 
previous situation. 

For a four-site ring system, if we only consider nearest 
neighbor couplings, we can easily get the analytical ex- 
pression of the residence time and the average trapping 
time. If a photon is absorbed at site 1, and the excita- 
tion is also trapped at site 1, the average trapping time 
is given by (t) = 4fc t ~ 1 . If a photon is absorbed at site 1, 
and the excitation is trapped at site 2 (or equivalently, 
site 4), the average trapping time is 



4 1 1 

{t} ~ h + f + r + h/2 



3(2r + k t ) 

8J 2 



(23) 



If a photon is absorbed at site 1, and the excitation is 



trapped at site 3, the average trapping time is 

4 T 3k t 
{t) = k t + Jl + ^- (24) 

By comparing the three equations for the four-site ring 
system, it is found that, the average trapping time of the 
first system is the smallest among the three scenarios, 
and that of the second is larger than that of the third. 

As indicated by the definition of (t), exciton popula- 
tion is proportional to the value of average trapping time. 
Thus we can also probe the dynamics of the system us- 
ing exciton population. It is interesting to note that the 
exciton population of site 3 grows larger than that of site 
2 and site 4 even it is far from site 1, where the exci- 
ton is initially located, the phenomenon also happens for 
the second and third scenarios. To explain this, one can 
consider the two distinct paths for the exciton produced 
at site 1, in clockwise and and counterclockwise direc- 
tions. They have the same phase when they meet at site 
3, which makes the exciton population greatly enhanced 
due to quantum interference. As a result, the exciton 
population of site 3 is larger than that of site 2 and site 
4 even though they are closer to site 1. 

Furthermore, we note a few differences among these 
scenarios : pi(t) is larger in the second scenario than 
others, therefore, it has the largest average trapping time, 
which leads to the lowest energy transfer efficiency, and 
P2{t) and Pi{t) are exactly the same in the first scenario 
due to positional symmetry. Pi(t) is slightly larger in 
situation three than in situation one, because the trap at 
site 2 destroys the symmetry. 



C. Energy Transfer in Single Ring Systems 

It is more interesting to look into the energy trans- 
fer efficiency of the ring system from the point of view 
of average trapping time, and as a starting point, simple 
single-ring systems with six, seven, eight and sixteen sites 
have been studied. The average trapping times at differ- 
ent trapping positions for the ring systems (with other 
parameters fixed to H = 1, kt = 2, T t = 20, J = 100) are 
given in Fig. [HI Site 1 is located at the origin (0°), and 
other sites are distributed in the ring uniformly. Taking 
a six-site system as an example, the zero-degree trapping 
corresponds to the case in which a photon is initially ab- 
sorbed at site one and the excitation is trapped at the 
same site. Similarly, (180°) trapping corresponds to the 
situation in which a photon is absorbed at site 1, and 
the excitation is trapped at site 4. It is interesting to 
find that the systems with an even number of sites have 
a comparatively smaller average trapping time when the 
exciton is trapped at the exactly opposite position (e.g., 
180°) of photon absorption. Since this phenomenon does 
not occur in systems with odd number of sites, it can be 
inferred from the simple ring systems that the construc- 
tive and destructive interference, closely related to the 
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FIG. 6: The average trapping time versus varying trapping 
positions for simple single-ring systems. From bottom up- 
wards, the square, circle, triangle, down triangle, left triangle 
and right triangle corresponds to three, four, six, seven, eight 
and sixteen-sites ring systems, respectively. For all systems, 
the parameters are chosen to % = 1, kt = 2, Ft = 20, and 
J = 100. 



symmetry of a system, plays a significant role in deter- 
mining energy transfer efficiency. 

For a ring system with even number of sites, it is pos- 
sible to have constructive interference at the site exactly 
opposite to the photon absorption site in the ring, where 
the exciton amplitudes are now in phase due to their 
symmetrical transfer path. In this scenario, construc- 
tive interference provides almost exactly the same aver- 
age trapping time compared to the case where exciton is 
created and trapped at the same site. However it should 
be noted that for systems with larger sizes, dephasing 
becomes increasingly important, resulting in a slight re- 
duction in the energy transfer efficiency. 

Now we turn to more realistic systems of photosynthe- 
sis, where capture of solar photons is often facilitated by 
light-harvesting antennae, and excitons, first created in 
the antennae, are transferred to reaction centers. 

Motivated by the results of simple ring systems, we 
examine further how the trapping position would affect 
excitonic energy transfer in a realistic ring system that 
includes dipolc-dipole interactions between all transition- 
dipole pairs that are not nearest neighbors. Fig. [TJa) 
displays the dependence of average trapping time on the 
trapping position, and the inset of Fig. [7f a) , a schematic 
of B850 ring of LH2, in which solid dots denote the a- 
apropotein, and the hollow dots, the /3-apropotein. Inter- 
estingly, a dip corresponding to a minimal average trap- 
ping time is found to appear at site 9, in addition to the 
one at site 1. In comparison to the uniform ring sys- 
tems shown in Fig. [6l there emerge zig-zag features in 
the trapping-time versus site- number plot in Fig. [7{a). 
Although the steepest dip can be simply explained with 
constructive interference, the appearance of multiple dips 
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FIG. 7: Average trapping time versus trapping site for the 
inset schematics of B850 ring of (a) LH2, and (b) LH1. Solid 
dots represent for the a-apropoteins and hollow dots the j3- 
apropoteins. The directions of the dipoles are tangential to 
the circle, with a-apropoteins pointing in clockwise direction, 
/3-apropoteins anti-clockwise. A photon is absorbed at la 
BChl, trapped at different BChls. For all calculations, the 
trapping rate kt is set to 1 and dephasing rate F is equal to 
5. The rest of parameters, i.e., J nm , are given in Eq. (0 and 
Eq. @. 



is directly related to the inclusion of long-range dipo- 
lar interactions, and therefore, multiple exction transfer 
pathways. Furthermore, it is also noted that the depen- 
dence of average trapping time on the site index no longer 
has the symmetries that exist for uniform rings. Further- 
more, the non-symmetric dependence of (t) on the trap- 
ping position makes it possible to have different clockwise 
and anti-clockwise paths, an effect that can be related to 
the dimerization of the LH2 ring. 

Dependence of average trapping time on trapping po- 
sition has also been examined for a single LH1 ring, and 
the results are shown in Fig.fTJb). The exciton Hamilto- 



nian of LH1 has a form similar to Eq. ([7]) with the total 
number of sites increased from 16 to 32. Chromophores 
are assumed to be evenly distributed with a distance of 
about 10 A, and other parameters for the LH1 system are 
taken from Ref. 2 -. Results similar to the LH2 case can 
be obtained from Fig. [7]Jb), and a dip appears around 
the middle point of the ring at site 17 with an obvi- 
ous loss of symmetry. Due to the increase in the total 
site number, which mitigates interference effects, a much 
smoother curve is obtained in Fig. [7f b) . For both the 
LH2 and LH1 rings, we have found that the long-range 
dipolar coupling reduces substantially the effect of quan- 
tum inference, in agreement with similar earlier findings 
in polaron dynamics simulation. 



D. Energy Transfer in Multiple Connected Ring 
Systems 

To expand our study of energy transfer efficiency to 
multiple-ring systems, two-ring toy models as depicted in 
Fig. [8] are examined first. Intra- ring nearest-neighbor in- 
teractions are assumed with coupling strength J, and site 
energies are set to be the same. The two rings are cou- 
pling through site 16 and site 17 with coupling strength 
V. Our findings are shown in Fig. [8j and an optimal 
coupling strength for efficient energy transfer is discov- 
ered at zero pure dephasing rate T = 0. This optimal 
strength found is close to the value of intra-ring coupling 
J, a result of system symmetry loss. If the trapping site 
is moved from site 32 to site 31, the optimal trapping 
time is found at a nonzero pure dephrasing rate T = 0.2, 
as shown in Fig. [8jc). It can be speculated that, to op- 
timize energy transfer efficiency, system symmetry im- 
perfections can be compensated by the introduction of 
pure dephasing or dissipation of system coherence by the 
environment. Such a finding is applicable for design of 
artificial photosynthetic systems in general. 

Next we will look at a realistic two-ring LH2 sys- 
tem with long-range dipolar interactions as illustrated 
in Fig. |9]Ja). Two rings interact via dipole-dipole cou- 
pling, according to Eq. ([5]). An initial state in which all 
sites in the first ring are evenly excited is selected, i.e., at 

t = 0, p nn = 1/16 for n = 1, 2, 16 in the first ring. We 

also assume that the excitation is evenly trapped in the 
second ring. It is found that the trapping time increases 
with the edge-to-edge interring distance r as expected, 
and for r > 20A (or with higher dephasing rates), the 
trapping time increases rapidly. It is known that quan- 
tum coherence may help enhance the efficiency of energy 
transfer 24 . At higher dephasing rates, with the destruc- 
tion of coherence, the excitation is transported classically, 
resulting in longer trapping times and lower efficiency. 

The findings on average trapping time in Fig.|H]suggest 
extremely distance-dependent nature of energy transfer. 
As indicated in Ref4A, the sites too close to each other in 
certain orientations will lead to non-fluorescent dimers, 
which can serve as excitation sinks that deplete excita- 
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FIG. 8: (a) Schematic of a two-ring toy model, (b) Contour 
map of (t) as a function of V and V, with trapping site locating 
at site 32 of the second ring as shown in (a), (c) Contour map 
of (t) as a function of T and V, with trapping site locating 
at site 31 of the second ring. The value of nearest neighbor 
coupling strength is set to J=l and the trapping rate kt = 1. 
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FIG. 9: (a) Schematic of a two-ring LH2 system, r is the 
distance between the edges of the rings, (b) Short distance 
(t) contour map as a function of T and r. (c) Long distance (t) 
contour map as a function of V and r. In the calculations of 
{t)s, the first ring on the left is evenly excited, and exicton is 
evenly trapped in the ring-2. Inter-ring interaction takes same 
form of the intra-ring dipole-dipole coupling. The trapping 
rate in the calculations is assumed to be a constant, with kt 
= 1. 



tions. 

One can also observe an optimal value of inter-ring 
distance of around 5A, corresponding to the highest 
strength of dipolar interaction between the nearest neigh- 
bor dipoles. This value is not close to the realistic inter- 
ring distance of LH2 complexes in a close packed mem- 
brane. Although it has been pointed out in Ref4i that 
a typical inter-chromophore center-to-center distances of 
neighboring molecules are consequently close ( 6 to 25 A). 
Thus, the global minimum obtained with the method in- 
troduced in this paper is only an ideal case such that two 
molecular aggregates are in contact with each other with- 



out any support of protein scaffold. Although the value 
only corresponds to the ideal case with highest inter-ring 
coupling strength, the existence of the local minimum is 
non-negligible. Such local minimum in average trapping 
time with non-zero dephasing rate appears in a continu- 
ous fashion throughout a large range of distances between 
the rings as shown in Fig. IHJb) and (c). The existence of 
this local minimum can be explained with the role of en- 
vironment related dephasing in increasing energy trans- 
fer efficiency for a system with asymmetric properties. 
The symmetry-breaking dipolar interaction in the LH2 
system leads to multiple non-nearest neighbor energy 
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transfer pathways, which can enhance the energy transfer 
within the ring greatly as shown in the polaron dynam- 
ics calculations, however it is more preferable to have 
more exciton populations that can be trapped at certain 
site. Thus it is preferable to have finite dephasing to 
increase the inter-ring transfer efficiency by maximizing 
the amount of exciton populations that can be trapped 
in the trapping site of the second ring through multiple 
pathway interference. This small amount of dephasing 
is not sufficient to destroy the quantum coherence in the 
system, which gives us a hint on how the unique arrange- 
ments and interactions between chromophores in the LH2 
complexes as well as the way that each complex interacts 
with others in optimizing the energy transfer in the pres- 
ence of an environment that is coupled to the system. 
More interestingly, in a recent work on Fenna-Matthews- 
Olson (FMO) complex^, even for this hnite-size, disor- 
dered molecular network can effectively preserve coher- 
ent excitation energy transfer against ambient dephasing, 
which lends a further support to what we have revealed in 
this section. Furthermore, the role of dipolar interaction 
in increasing exciton transfer efficiency through introduc- 
ing of multiple pathway for inter-ring exciton transfer is 
evident, which further supports the previous findings. 



for the coupled rings, an optimal coupling strength is ob- 
tained between the hypothetic two-ring systems with zero 
dephasing rate, while in a two-LH2 system, owing to the 
intrinsic symmetry breaking, global minimum and local 
minimum of average trapping time is always found with 
non-zero depahsing rate. In both the cases, the efficiency 
drops with further increase of dephasing rate, due to the 
loss of quantum coherence. The findings in this paper re- 
garding exciton transfer processes in ring systems might 
assist us to design more efficient artificial light harvesting 
systems. 
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Appendix A: Precision of the Time-dependent 
Davydov Ansatze 



IV. CONCLUSIONS 

In this paper, two distinct approaches have been ap- 
plied to study the excitation energy transfer in the ring 
systems. From both methods, non-nearest neighbor dipo- 
lar interaction is found to be helpful in increasing intra- 
or inter-ring exciton transfer efficiency as a result of mul- 
tiple pathways. 

Using the Haken-Strobl model, the energy transfer in 
the ring systems with coupling to environment is stud- 
ied in details. The energy trapping efficiency is found 
to be dependent on the trapping position in both the 
hypothetic and the realistic ring systems. For the hy- 
pothetic system, with the trap positioned diametrically 
opposite to the site of the initial excitation (which is 
possible only in systems with even number of sites) , we 
observe a sudden drop of the trapping time due to the 
constructive interference between the clockwise and anti- 
clockwise paths leading to the trap. The situations in a 
LH2 (B850) ring and LH1 ring are similar to the hypo- 
thetic systems, but owing to the broken symmetry caused 
by the dimerization of BChls and dipolar couplings, the 
drop of the trapping time is not significant. In other 
words, the ring becomes more homogeneous in terms of 
energy transfer, which further supports the findings of 
polaron dynamics studies of one-ring system. 

The polaron dynamics of the coupled rings revealed 
interesting role of dipolar interaction in increasing en- 
ergy transfer efficiency by introducing multiple transfer 
paths between the rings, and the time scale of inter-ring 
exciton transfer is in good agreement with previous stud- 
ies. Moreover, when the Haken-Strobl model is applied 



For a trial wave function that does not strictly 

obey the time-dependent Schrodinger equation, the de- 
viation vector \6(t)) can be defined as 



|«f(t))=t^|*(t))-£|*(t)), 
and the deviation amplitude A(£) is defined as 



A(t) EE vWlf. 



(Al) 



(A2) 



For the Holstein Hamiltonian H defined in Eqs. (JTJ- 
0, one can derive the explicit expression of (5(t)\5(t)) 
for the Di Ansatz: 



\iu n {t) + T n {t) + a n (t)R n (t) 

n 

9q 

a n (t)jta q 



(A3) 



with 



R n (t) = Re^ [iX ntq {t)+^Lo q e-^ n ~uj q X^ q {t)]X^ q {t), 



N 

(A4) 

where T n (t), Q, n ,q{t) and Ajj(t) are three expressions 
which have no item of a n (t) or X n ,q{t)- 
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For the D2 Ansatz, Eq. (|A3I) can be further derived to 

(w*)|w*)) 

= J2\ia n (t)+T n (t) + a n (t)R n {t) 



+ Y,\0 q (t) + jl=u q J2\<Xn(t)\ 2 e 

q n 



1 V - ^ 9 2 



i-|Ekmi 



2 £ -iqn 



tqn -LO q P q (t) 

(A5) 



And for the D Ansatz, Eq. (|A3|) can be further derived 
to 



v *6(*)|*d(*) 
= E \ia n (t) + T n {t) + a n (t)R n (t) 



+ \iXo(t) + -7=^0 ~ w A (i) 
v TV 



l"^(t)| a 



E EK«w| - i»«(*)i s 
a dW-EE|"^w 



KWI 2 



i-|^(*)l s 



(A6) 



where £ g (t), 0n ig (i), j/ g (t) and d q {t) are four expressions 
which have no item of a n (t), $ q (t), j q (t) or Xo(t). 

Substituting Eq. (JTrj into Eqs. one obtains 

the expressions for the system energies by the Davydov 
Ansatze in the Holstein model: 

E„(t) = (*D(t)|flBx|*D(t)) 

- -2JReE<(*)Sn,n+i {t)a n+1 (i), (A7) 



£ ph (t) = (* D (i)|i? ph |* D (t)) 

= e[kwi 2 E^iwoi 2 ], (as) 



and 



£ cx _ ph (t) = (* D (t)|fl-„_ph|* D (t)) 

2 E [|«„W| 2 ReE5 9 ^A„, g (t) e ^« 



iV 



where S n>m (t) is the Debye- Waller factor. 



(A9) 



Note that since the unit of A(t) is that of the energy 
by comparing A(t) with the main component of the 
system energies such as E p h(t) and £ , ex _ p h(t), one can 
observe whether the deviation of an Ansatz from obeying 
the Schrddinger equation is negligible or not, in the 
concerned case. From this perspective, the comparison 
between A(i) and energy components of the system 
provides a good reference for the validity of an Ansatz. 



Appendix B: Detail derivation of polaron dynamics 
in multiple-ring system 

The system energies of the multiple-ring system can 
be obtained as follows by apply Davydov Ansatze to the 
modified Holstein Hamiltonian given in Eq. ([5]), © and 
CD): 

E cx (t) --EE J ™ <*W«m (*)^m (*) (Bl) 



Bp h («) = EEK(*)| 2 E^I A n 9 (i)| 2 (B2) 

r n q 
* r n q 

(B3) 

with Debye- Waller factor S^]^ 2 (t) given by: 

S™»(t) = exp{E[A^(t)A^(t)5 rir2 

- \\K\{t)? - ^|A^(«)| 2 ]}, (B4) 

and the index r% and r 2 all runs over 1 to N T - lng . 

Since Dl is the used for all calculations in this paper, 
here we only give the expression of (8(t)\6(t)) for Dl, 
which can be derived as follows: 

<<W*)I<W*)) 

r n 

+ E E E +n »*(*)i a+A ^*) 

r o£(t)#0 q 

-EE Ek^-EE^i 2 - (^5) 



Minimization of the first two terms in Eq. (|B5|) lead 
to the equations of motions for the time-dependent vari- 
ational parameters ct r n {t) and \ r nq (t) as: 



l (t)=i[T^t) + a r n (t)R r n (t)} ) 



(B6) 



and 



fl r (t) a r 

= ^~TTiY + J h< e ' iqn ~ < x nM- ( B7 ) 



14 



Each collected term in Eqs. (IB6|) and (IB7|) is given as: 

< g (t) = d r m (t) + ^>~ iqn - <KM (B8) 



=EEir<(t)C(«)[AS ? (*)^-A; i ,(t)], 



and 



(Bll) 



i£(f) = ReJ2K q (t) + -^<e-^"]A-(t), (B9) 

8 * 



(BIO) 



a dW = EEtE E^jsr^w^w^wi 



fl n rim' r 2 m 



(B12) 
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